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ABSTRACT 

Extensive photometric and spectroscopic surveys of the Andromeda galaxy (M31) have discovered tidal debris 
features throughout M31's stellar halo. We present stellar kinematics and metallicities in fields with identified sub- 
structure from our on-going SPLASH survey of M31 red giant branch stars with the DEIMOS spectrograph on the 
Keck II 10-m telescope. Radial velocity criteria are used to isolate members of the kinematically-cold substructures. 
The substructures are shown to be metal-rich relative to the rest of the dynamically hot stellar population in the fields 
in which they are found. We calculate the mean metallicity and average surface brightness of the various kinematical 
components in each field, and show that, on average, higher surface brightness features tend to be more metal-rich 
than lower surface brightness features. Simulations of stellar halo formation via accretion in a cosmological context 
are used to illustrate that the observed trend can be explained as a natural consequence of the observed dwarf galaxy 
mass-metallicity relation. A significant spread in metallicity at a given surface brightness is seen in the data; we show 
that this is due to time effects, namely the variation in the time since accretion of the tidal streams' progenitor onto 
the host halo. We show that in this theoretical framework a relationship between the alpha-enhancement and surface 
brightness of tidal streams is expected, which arises from the varying times of accretion of the progenitor satellites 
onto the host halo. Thus, measurements of the alpha-enrichment, metallicity, and surface brightness of tidal debris 
can be used to reconstruct the luminosity and time of accretion onto the host halo of the progenitors of tidal streams. 
Subject headings: galaxies: substructure — galaxies: halo — galaxies: individual (M31) 



1. INTRODUCTION 

In a hierarchical universe, a large galaxy will undergo 
many mergers during its lifetime. The amount and prop- 
erties of accreted material observed in present-day galax- 
ies can be used to test hierarchical formation scenarios. 
Stellar halos provide an ideal environment for investigat- 
ing the detailed merger history of an individual galaxy. 
The sparse stellar populations of a galaxy's halo com- 
bined with long dynamical times make it possible for 
tidal debris features to remain identifiable in phase-space 
for billions of years. In contrast, it is difficult to sepa- 
rate material formed in situ from accreted material in a 
galaxy's relatively dense disk and bulge. 

Recent sophisticated numerical and semianalytic sim- 
ulations of stellar halo formation have made great 
strides in characterizing the properties of stellar ha- 
los that are built up through the tidal st r ipping of 
accre t ed sa tellites (e.g., i Johnston et al.l [T996I: iJohnstonI 
lOM iHelm i & White 199 91 iHelmi fc de Zeeuwl 120001: 
Bullock et al...2001i ). iBuUock fc JohnstonI (|2005[ l (here- 
after BJ05) used a combination of iV-body and semi- 
analytic techniques to simulate the formation via ac- 
cretion of 11 MW-analog stellar halos in a cosmologi- 
cal context. These high resolution simulations allow de- 
tailed analysis of the expected physical properties of stel- 
lar halos composed primarily of tidal debris and their 
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associated satellite systems ([Font et al.l l2006al |bl. 120081 : 
IJohnston eralll2008[ ). 

The stellar halos of the Milky Way (MW) and An- 
dromeda (M31) galaxies both show ample evidence 
of recent and on-going accretion events. Among 
the most prominen t of these substructures are the 
Sagittarius stream (Ibata et al.l 119941 : iMajewski et al.l 

t03; Ncwbcrg ct al. 2003) and the M onoceros stream 
annv et al.ll2003t Rocha-Pinto ct al. 20 031) in the MW , 
and the giant southern stream (CSS; Ibata et al.ll200ll) 
in M31, which likely po llutes much of M31 ' s spheroid in- 
terior to ~ 20-30 kpc (iFardal et aLl l2007l: iGilbert et all 
I2007t iRichardson et al.ll2008f) . 

In addition to these prominent stellar streams, re- 
cent large observational surveys are increasing the num- 
ber of known tidal debris features by pushing to lower 
surface brightness limits. Analysis of the Sloan Dig- 
ital Sky Survey has revealed many fa int tidal debris 
features in the MW 's stellar hal o (e.g., [G rillmair 2006|; 
iGrillmair fc Dionatos 2006: Be lokurov et a l. 2006, 2007), 
while large-scale photometric studies of M31 arc reveal- 
ing similarly large numbers of tidal debris features in 
M31's stellar halo (|Ibata et all |2001L 1200?!) . Spectro- 
scopic survey s of re d giant branc h (RGB ) stars in M31 
(llrwin et all 120051 : iKalirai et all l2006at IGilbert et all 
l2006( l are enabling the discovery of tidal debris fea- 
tures thrwoghthe kinematic s of M31's stellar p opulations 
(e.g.JIbata eTai . 2005; Kalir ai_eiall20 06b: Gi lbert et all 



120071 IChapman et al.li200a : IGilbert et al.ll2009l ). 



A classic signature of substructure is an abnormally 
high surface brightness relative to the value expected 
from the observed surface brightness profile of the 
stellar halo. Due to the long dynamical times in 
the sparse outer regions of galaxy halos, tidal debris 
from relatively recent mergers (within the last several 
Gyr) will not have had time to become spatially well- 
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mixed. The kinematical distribution of a well-mixed 
stellar halo population is expected to be broad (dy- 
namically hot) out to large distances from the cen- 
ter of the galaxy. In the MW and M31, the mea- 
sured line-of-sight velocity dispersion i s ay ^ 120 to 
150 km s"^ near the galax y 's center (iBattaelia et al.l 
l2005t iChapman etlll [200l [Gilbert et all l2007fl and 
decreases to ^ 1 00 km s~^ at R ~ 60 kpc 
(jBattagUa 'et all 120051 : IChapman et all 120061 ). In con- 
trast, tidal debris from a single accretion event should 
be dynamically cold, since the local velocity dispersion of 
the debris should decrease as it mixes along the progen- 
itor satellite's orbit: as the stars spread apart, decreas- 
ing their spatial density, their density in velocity-space 
must increase in order to conserve t heir phase-space 
density and satisfy Liouville's theorem (jHelmi fc Whitd 
|1999| ). Hence, kinematical substructure in a survey is 
another signature of hierarchical structure formation. A 
difference in chemical abundance between stars in the 
kinematically hot and cold populations or between two 
closely spaced fields is also evidence of substructure, as 
the progenitors of recent accretion events are likely to be 
chemically distinct from the ensemble of satellites that 
were accreted at early times. 

We have utilized the kinematical and chemical signa- 
tures of substructure to iden tify multiple, ki nematically 
cold tidal debris features (|Guh athakurta c t al.l l2006t 
iKalirai et al1l2006bl : [Gilbert et al..,2007..2009.) during the 
course of the Spectroscopic and Photometric Landscape 
of Andromeda's Stellar Halo (SPLASH) Survey which 
includes extensive Keck/DEIMOS^ spect roscopic survey 
of M31's stellar ha l o fe.g.. [Guh athaku rta et all 120051 : 
iKaUrai et all l2006at iGilbert et al.. ,2006 . l2007f ) Kine- 
matical signatures of substructure enable the detection 
of faint and diffuse features which are not easily iden- 
tified in sta rcount maps. For example, a shelf feature 
([Gilbert et all l2007h along M31's southeast minor axis 
(JSproj ~ 12-18 kpc) went unnoticed on star count 
maps (in fact, the region was specifically targeted for 
deep HST photometry based on its apparent smoothness; 
I Brown et a l. 2003). The observed distribution of stars in 
line-of-sight velocity vs. position space and its striking 
similarity to model predictions led to the discovery of 
the SE shelf. 

We can use the kinematically cold tidal debris features 
found in our survey to study general trends of tidal de- 
bris in stellar halos. The recent advances in simulations 
of stellar halo formation have made it possible to under- 
take detailed comparisons o f the properties of observed 
and simulated substructure (jBell et al.l 120081 : iFont et all 
l2008l ) . In this contribution, we present a relationship be- 
tween the surface brightness and metallicity of tidal de- 
bris features observed in our M31 data set, and provide 
a physical interpretation of the observed trend through 
comparisons of the observations with simulated stellar 
streams in the BJ05 simulations. 

The paper is organized as follows. In §[2] we briefly 
discuss the observations and reduction techniques and 

^ Data presented herein were obtained at the W. M. Keck Ob- 
servatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California and 
the National Aeronautics and Space Administration. The Obser- 
vatory was made possible by the generous financial support of the 
W. M. Keck Foundation. 



present the fields in our M31 survey in which substruc- 
ture has been identified. We present the kinematical and 
chemical abundance distributions of the stars in these 
fields as well as the surface brightness and average metal- 
licity of the observed kinematical substructures. The 
methods used in the BJ05 simulations are outlined in 
§131 and the trends of surface brightness with chemical 
abundance of stellar streams in the simulated BJ05 stel- 
lar halos are presented. We compare the observations 
with the simulations and discuss the resulting physical 
implications in §|4l Finally, we summarize our results in 

m 

2. DATA 

2.1. Photometric and Spectroscopic Observations 

The photometric and spectroscopic observations and 
data reductions employed in our SPLASH M31 sur- 
vey h ave been discussed in d e tail in previous pa- 
pers i Guhathakurta et all 120061: iKalirai et all l2006bl : 
IGilber t ct al. 2006, 200'5^and will be only briefly out- 
lined here. The observations are organized into "fields" 
in M31's stellar halo based on the spatial proximity of 
the spectroscopic slitmasks (Table |T]) . A single photo- 
metric pointing is large enough to support multiple, non- 
overlapping DEIMOS slitmask observations; the spectra 
from slitmasks based on a single photometric pointing 
are in general analyzed together and considered a field. 
The exception to this are fields f207 and H13s, which 
are from a single photometric pointing but are treated 
as separate fields due to the fact that they are located 
at different positions along M31's GSS; the debris from 
the GSS is centered at a different mean velocity in each 
of these fields. 

Photometry and astrometry for the SE shelf, f207, and 
H13s fields were derived from images in the g' and i' 
bands taken with the MegaCam instrument on the 3.6- 
m Canada-France-Hawaii Telescope (CFHT).^ The SEx- 
tractor program was used for object detection, photom- 
etry, and morphological classification (jBertin fc Arnouti 
Il996f) and observations of Landolt photometric standard 
stars were used to transform from instrumental g' and 
i' magnitudes to Joh nson-Cousins V and / magnitudes 
(jKalirai et al.ll2006bD . 

Photometry and as trometry fo r the a 3, al3, and m4 
fields were derived bv lOstheimeil ()2003l ) from images in 
the Washington M and T2 bands and the intermediate- 
width DD051 band taken with the Mosaic camera on the 
4-m Kitt Peak National Observatory (KPNO) telescope.^ 
The DD051 filter includes the surface-gravity sensitive 
Mg& and MgH stellar absorption features. The combina- 
tion of these three filters allows photometric selection of 
stars that are likely to be M31 red giant branch (RGB) 
stars rat her than MW dwarf s tars along the line-of-sight 
to M31 (|Maiewski et al.ll2005D. The photorn etric trans- 
formation relations in iMaiewski "eTall ((2OOOI) were used 

MegaPrime/MegaCam is a joint project of CFHT and 
CEA/DAPNIA, at the Canada- France-Hawaii Telescope which is 
operated by the National Research Council of Canada, the Insti- 
tut National des Science de I'Univers of the Centre National de la 
Recherche Scientifique of France, and the University of Hawaii. 

^ Kitt Peak National Observatory of the National Optical As- 
tronomy Observatory is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 
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TABLE 1 

M31 Spectroscopic Fields with Identified Substructure. 



Field 


# Spectroscopic 


Projected # Stars^ 


# of M31 Reference 




Masks 


Radius 


Stars'^ 






(kpc) 





SE shelf 


4 


12- 


-18 


510 


428 


Gi 


bert et a 




f207 


1 


16- 


-19 


57 


49 


Gi 


bert et a 





H13s 

a3 

al3 
m4 



20-23 

31-34 

55-63 
53-60 



212 



109 
154 



178 

69 

44 
53 



Kalirai et al 
Gilbert ct al 
Guhatliakurta ct al, 
Gilbert ct al. (2009) 
Glbcrt ct al. (2009) 
Gilbert et al. (2009) 




^ The total number of unique stellar spectra for which we recovered velocities.^ The number of M31 
RGB stars is defined as the number of stars that are identified as secure and marginal M31 RGB stars 
by the [Gilbert et al.l l l2006f) diagnostic method. 



to convert the M and T2 magnitudes to Johnson-Cousins 
V and / magnitudes. 

Objects were selected for spectroscopy based on their 
derived / magnitudes and morphological properties. In 
the outer fields (a3, al3, and m4), objects were further 
prioritized by their probability of being an M31 RGB 
star based on the M, T2, and DD051 photometry; this 
is vital for increasing the spectroscopic efficiency of our 
survey in the sparse outer regions of M31's stellar halo. 
Spectroscopic slitmasks were observed with the DEIMOS 
instrument on the Keck II telescope for 1 hr each, using 
the 1200 line mm^^ grating (which has a dispersion of 
0.33 A pixer 1) centered at 7800 A. This led to an average 
spectral coverage of 6450-9150 A. 

Modified versions of the spec2d and specld soft- 
ware* developed by the DEEP2 team at the Univer- 
sity of California, Berkeley were u sed to reduce and an- 
alyze the spectrosc opic data (see ISimon fc Gehal [20071: 
iGilbert et aLll2007t ). Reduced one-dimensional spectra 
are cross-correlated with a library of template stellar 
spectra to determine the redshift of the object. Each 
spectrum is visually inspected and assigned a quality 
code based on the number and quality of the absorp- 
tion lines. A spectrum must have at least two identi- 
fiable spectral features to b e considered to have a se- 
cure redshift meas urement (|Guhathakurta et al.l l2006t 
IGilbert et al.lf2007t ) . A heliocentric correction and a cor- 
rection for possible miscentering of the star within the 
slit (based on the position of the telluric A-band relative 
to atmospheric emission lines in the obse rved spectrum) 
is applied to the mea sured velocities (jSimon fc Gehal 
I2007t ISohn et all 120071) . The median velocity error is 
~ 6 km s~^, based on error estimates from the cross- 
correlation routine and repeat measurements of individ- 
ual stars. 

2.2. Selection of M31 RGB Stars 

The line-of-sight velocity distributions of foreground 
MW dwarf stars and M31 RGB halo stars overlap, mak- 
ing it nontrivial to identify individual stars as M31 red 
giants or MW dwarfs. We use the diagnostic method 
described in [Gilbert et all (|2006l) to isolate M31 RGB 
stars from foreground MW dwarf stars in our spectro- 

^ http://astron.berkelev.edu/~cooper/deep/spec2d/primer.ht ml. 
|http://astron.berkeley.edu/~cooper/deep/specld/primer.html| 



scopic sample. Empirical probability distribution func- 
tions in five photometric and spectroscopic diagnostics 
are used to determine the likelihood an individual star is 
an M31 red giant. These diagnostics include (1) line-of- 
sight velocity, (2) photometry in the M, T2, and DD051 
bands (when available), (3) the EW of the Nai doublet at 
8190 A versus the {V — I)q color of the star, (4) position 
of the star in an [/q, {V — /)o] color magnitude diagram, 
and (5) photometric versus spectroscopic metallicity es- 
timates. The M31 RGB sample used in this analysis 
includes all stars that are more probable to be M31 red 
giants than MW dwarfs. 

2.3. Metallicity Estimates 

Metallicity estimates ([Fe/H]) are derived by compar- 
ing a star's position in the [Iq, {V — T)o] color mag- 
nitude diagram with theoretical isoch rones adjusted 
to the distance of M 31 (783 kpc; Stanek fc GarnavichI 
119981 : |Hollaiid|[i99 l). Interpolation within the grid of 
12 Gyr, [a/Fe]=0, IVandenBerg et al.l (|2006( l isochrones 
yields an [Fe/H] me asurement for each M31 RGB star 
(|Kalirai et al.ll2006aD . 

Spectroscopic estimates of the metallicity of an in- 
dividual star can also be obtained using the equiva- 
lent width of the Call t ri plet an d the calibration re- 
lations of iRutledge et"aI1 (|1997al |bh. On average, the 
two estimates are in good ag reement for the M31 sur- 
vey data (Kali rai et al.|[20b6aD . although there is signifi- 
cantly larger scatter in the spectroscopic metallicity esti- 
mates due to l a,rge measurement errors in the Ca 11 EW 
measurements (iKalirai et al.l[2006al : [Gilbert et al.l[2006h . 
This is due primarily to the relatively low average signal- 
to-noise ratio of M31 RGB spectra in our spectroscopic 
survey (typical values of S/N ~ 10 pixel"^). There is 
smaller scatter in the spectroscopic metallicity estimates 
and very good agreement with photometric metallicity 
estimates for our highest S/N M31 RGB spectra. Due to 
the large measurement uncertainties in the spectroscopic 
metallicity estimates for individual stars, we choose to 
use photometric metallicities for this analysis. 

The use of photometric metallicities introduces system- 
atic biases due to the isochrone set used and to the adop- 
tion of a single age and alpha-enrichment for the stellar 
populations. Estimates of [Fc/H] based on different sets 
of isoc hrones vary by less than 0.15 dex ([Kalirai et al.l 
[20063) • Although the alpha-enrichment of M31's stel- 
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lar halo has yet to be observationally constrained, deep 
HST/ACS observations of fields in M31's stellar halo 
which reach the main sequence turnoff fe.g. jBrown et al] 
l2006al |2007[ ) have shown that the stellar population of 
M31's halo spans a wide range in age. In the SE shelf 
and H13s fields, which are contaminated by GSS debris, 
the stellar population ranges from 6-13 Gyr in age. The 
younger stars are found to be the most met al-rich, while 
the old population is relatively mctal-poor ([Brown et al.l 
I2006bllaf ). In an HST/ACS M31 halo field at 21 kpc, a 
spread in ages of ~ 8-13 Gyr is found. Once again, 
the younger stars are the most metal-rich while the old 
stellar population is relatively metal-poor. We therefore 
expect our metallicity estimates to be affected mostly at 
the metal-rich end by the assumption of a single, old stel- 
lar population. Varying the age from 12 to 6 Gyr results 
in an « -f0.15 dex shift in the [Fe/H] estimates derived 
from isochrone fitti ng for metal-rich ([Fe/II]> —1) stars 
(|Gilbert et al.l[2009 l 

2.4. Identification of Substructure: Kinematics and 
Metallicity 

Figure [T] presents line-of-sight velocity distributions 
and [Fe/H] versus line-of-sight velocity for the stars in 
each of six fields in which our survey of M31's stellar 
halo has discovered kinematical substructure. The distri- 
butions of CMD-based [Fe/H] values (§[231) in Figure [T] 
show the presence of a relatively metal-poor ([Fe/H]< 
— 1) population with a broad velocity distribution in 
each field, while the kincmatically cold substructure 
tends to be relatively metal-rich ([Fe/H]> —1). Detailed 
analyses of the kinematics, chemical properties, spatial 
distributions, and plausible origins of the various sub- 
structures are prese nted el sewhere [SE shelf region i n 
[Gilbert et all (|2007|): a3 in iGuhathakurta et all (j2006D : 
H13s in lKalirai et all (l2006bl): fields 1207. H13s. a3. al3, 
and m4 in [Gilbert et all (|2009f) ] . 

The primary signature used to identify substructure 
in a field is the presence of a kinematically cold pop- 
ulation; the line-of-sight velocity distributions of these 
fields are inconsistent with being drawn from a sin- 
gle, broad Gaussian centered at M31's s ystemic veloc- 
ity (( t))^- 300 km s-\ a-u = 129 km s~^: [Gilbert et all 
l2007t ). Fields al3 and m4 are at larger projected 
radial distances from the center of M31 (i ^proi ~ 
60 kpc) than the study of [Gilbert erall ([2007[ ). Their 
line-of-sight velocity distributions are also inconsistent 
with being drawn from a Gaussian with parameters 
300 km s^^, cr^u = 100 km s~^, the velocity disper- 
si on implied for th e dista nce of these fields by the study 
of IChapman et al.l ([20061 ). 

The velocity distribution in each field has been fit by 
a combination of Gaussians using a maximum-likelihood 
technique (red curves in Figured]). In general, the num- 
ber of Gaussians used in the multi-Gaussian fits is the 
minimum required to produce a maximum-likelihood fit 
that is consistent with the data, although additional ev- 
idence of substructure in a field is also taken into ac- 
count in determining the number of kinematical compo- 
nents in a given spectroscopic field ([Gilbert et al.|[200"7l 
l2009t ). The Gaussian components of the underlying 
kinematically hot population were held fixed {{v^p^ = 
-300 km s-\ alP^^ 129 km s'^) while the prop- 
erties of the Gaussians representing the kinematically 
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Fig. 1. — Line-of-sight velocity distributions and [Fe/H] versus 
line-of-sight velocity for the fields discussed in this paper. Shaded 
velocity histograms and open circles represent stars classified as 
M31 red giants by the [Gilbert et al.l (2006) diagnostic method, 
while crosses represent stars classified as MW dwarfs. Open ve- 
locity histograms show the velocity distribution of all the stars 
(M31 RGB and MW dwarf) in our fields. The solid curves show 
the best-fit velocity distributions in each field (§[2]4]l. In general, 
the observed kinematically cold populations arc more metal-rich 
than the rest of the stellar population. 



TABLE 2 

Fraction of M31 RGB St.-^rs in Kinematically Cold Components. 



Field 


( km s"-"-) 


AfKCc/A^tot'' 


SE shelf 


-285 


41+"-^" 


f207 


-524 


Q O1+0.08 
'-'■■'^-0.07 


f207 


-426 


" 
"■•^^-0.10 


H13s 


-490 


^■^'^-0.06 


HlSs 


-388 




a3 


-441 


0.59t°:l^ 


al3 


-301 


'^-0.21 


m4 


-355 


45+0 " 


m4 


-255 


16+°-l'' 



cold populations (mean velocity and velocity dispersion) 
and the fraction of stars in each cold population were 
allowed to vary. The SE shelf, f207, H13s, a3, and 
al3 fields all contain substr ucture related to the GSS 
in M31 tKalirai et al. 2006bl : [Guhathakurta et al.[ [20061 : 
[Gilbert et all 120071 l20'09h . Fields 1207 and H13s both 
have a second kinematically cold component which may 
be unrelated to the GSS. The cold components with 
mean velocities of —527 km s~^ and —490 km in 
f207 and H13s, respectively, are those associated with 
the GSS ([Gilbert et al.|[2009l ). Although there are some 
lines of evidence that suggest the substructure in m4 
and the secondary substructures in fi elds f207 and HI 3s 
ma y also b e associated with the GSS (Far dal et al.ll2008l : 
Gilbert et al. 2009), the origins of these features are cur- 
rently unknown. 
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Fig. 2. — Surface brightness versus average metallicity for each of 
the kinematical components in Figure [T] kinematically cold com- 
ponents likely associated with the GSS (solid circles), kinematically 
cold components which are not known to be associated with the 
GSS (solid triangles), and the kinematically hot components (open 
squares). The error bars represent the Poisson error for the sur- 
face brightness estimates and the error in the mean [Fe/H] value. 
The mean and root- mean-square (rms) [Fe/H] and fiY of the lower 
surface brightness {fiy > 29 mag arcsec"-^) and higher surface 
brightness (fiy < 29 mag arcsec"^) kinematically cold compo- 
nents are also shown (open hexagons and error bars). Higher sur- 
face brightness features are on average more metal-rich than lower 
surface brightness features, and kinematically cold substructures 
are on average more metal-rich than the kinematically hot popu- 
lations. Lower surface brightness features display a larger spread 
in metallicity than higher surface brightness features. The dotted 
line denotes the faintest identifiable tidal streams to date, given 
the number of stars in our DEIMOS survey. 



2.5. Characteristics of Observed Substructure: Surface 
Brightness versus Metallicity 

Figure[2]d.isplays surface brightness versus mean metal- 
licity for the 1VI31 stellar halo data presented in Fig- 
ure [T] Each kinematical component identified in the 6 
fields shown in Figure [1] is represented: the kinemati- 
cally cold components related to the GSS (soHd circles), 
the kinematically cold components that are not known 
to be associated with the GSS (solid triangles), and the 
kinematically hot population (open squares). 

Table [2] lists the maximum-likelihood fraction of 1VI31 
RGB stars in each of the kinematically cold components 
(from the fits shown in Fig. [T|). These fractions were 
used to determine the number of M31 RGB stars in each 
component. The ratio of M31 RGB to MW dwarf stars, 
multiplied by the expected number of JVIW dwarf stars 
along the line of sight to each fiel d (based on the B e- 
sancon Galactic population model; iRobin et al.ll2003t ) is 
used to estimate the surface brightness of each compo- 
nent, as described in iGuhathakurt a et al.l ()2005[ ). Cor- 
rection factors that account for the rate of successful 
velocity measurements as a function of magnitude and 
varying 1VI31 RGB target efficiency due to photomet- 
ric pre-selection of stars likely to be red giants are ap- 
plied to each field (Gilbert et al., in preparation). The 
spectroscopic surface brightness estimates are normal- 
ized using photometric sur f ace br ightness estimates from 
iPritchet fc van den Berghl (|1994[ ). Error bars represent 



Poisson error estimates based on the number of stars in 
each component. 

The [Fe/H] for the kinematically cold components is 
the mean [Fe/H] of stars within ±2av of the mean ve- 
locity of the cold component (jwstar- < 2(t^^^), 
while the [Fe/H] for the kinematically hot components 
is the mean [Fe/H] of stars with velocities greater than 
mean velocity of the cold components 
([wstar- (w)^^^ I > 2 cr^^^; in fields with two cold compo- 
nents, stars must meet this criterion for both). The num- 
bers of stars which are thus used for the [Fe/H] measure- 
ment for each component are consistent with the number 
of stars expected to be in each component based on the 
maximum-likelihood fits shown in Figure [T] Error bars 
represent the error in the mean [Fe/H] value. 

Since stars cannot be identified individually as belong- 
ing to a particular kinematical component, these veloc- 
ity cuts are chosen to maximize the number of stars in- 
volved in the [Fe/H] measurement of the kinematically 
cold component, while minimizing contamination in the 
mean [Fe/H] measurement from other kinematical com- 
ponents. Even so, the metallicity estimate of the sub- 
structure will be biased by contamination from the kine- 
matically hot component. This effect will be greatest for 
the least dominant substructures, and will in general bias 
the mean metallicity of the kinematically cold substruc- 
tures to lower values. In the SE shelf fields, where we 
have a large statistical sample, accounting for the con- 
tamination by the kinematically hot component in the 
measurement of the mean metallicity of the kinemati- 
cally cold comp onent changes the mean metallicity by 
only -1-0.03 dex (jGilbert et al.ll2007l ): however, it must be 
noted that the [Fe/H] distribution of the kinematically 
hot population in the SE shelf fields is very similar to the 
[Fe/H] distribution of the kinematically cold population. 

The assumption of a single, old stellar population likely 
provides a larger source of systematic error in the photo- 
metric (ClVID-based) metallicity estimates. As discussed 
in ii l2.3i this assumption is known to be wrong for at 
least some of our fields and is expected to affect our 
metallicity estimates mostly on the metal-rich end. In- 
deed, a comparison of the mean Call-based spectroscopic 
metallicity estimates (measured in bins of photometric 
metallicity) of stars associated with the kinematically hot 
and cold populations shows that the kinematically cold 
populations do, on average, have a higher spectroscopic 
abundance than the kinematically hot population at a 
given photometric abundance. This indicates that the 
substructure in JVISl's stellar halo is on average younger 
and/or more enriched in [a/Fe] than JVISl's kinemati- 
cally smooth stellar halo. Using 6 Gyr instead of 12 Gyr 
isochrones changes the [Fe/H] estimates by ^ -t-0.15 for 
more metal-rich stars ([Fe/H]> —1.0) and by as much 
as ~ -1-0.3 for metal-poor stars ([Fe/H]< —1.5). Thus, 
the true metallicity spread is likely greater than Figure [2] 
suggests. 

The data in Figure [2] show a clear correlation be- 
tween surface brightness and metallicity: higher surface 
brightness features are more metal-rich than lower sur- 
face brightness features. The kinematically cold sub- 
structures are on average more metal-rich than the un- 
derlying kinematically hot components, as was implied in 
Figure [T] Figure [5] also shows evidence for a larger spread 
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Fig. 3. — The same as Figure [2] but for entire fields from our 
M31 spectroscopic survey. Fields are divided into three categories: 
fields with (solid circles; Fig. [TJ and without (open triangles) kinc- 
matically identified substructure, and fields with small numbers 
of M31 RGB stars (< 10), for which we cannot determine if sub- 
structure is present. The large symbols show the mean and rms 
(errorbars) [Fe/H] and /ly of fields with and without identified 
substructure. Fields with identified substructure are on average 
more metal-rich than fields without identified substructure. The 
dashed line represents the current surface brightness limit achieved 
in our M31 spectroscopic survey, while the dotted line denotes the 
faintest tidal stream identified in our survey (Fig. [2j. 



in metallicities among lower surface brightness features 
{nv > 29 mag arcsec"^) than higher surface brightness 
features (/xy < 29 mag arsec"^). The root-mean-square 
(rms) deviation of average [Fe/H] values among higher 
surface brightness, kinematically cold components is only 
0.08 dex, while the rms deviation of average [Fe/H] val- 
ues among lower surface brightness, kinematically cold 
components is 0.32 dex. 

Figure[3]displays surface brightness versus mean metal- 
licity for fields from our M31 halo spectroscopic sur- 
vey with kinematically-identified substructure (those 
presented in Figure (TJ solid circles), fields without 
kinematically-identified substructure (open triangles), 
and fields for which we cannot determine if substructure 
is present due to small numbers of M31 RGB stars (< 10 
M31 stars per field; open squares). The fields with small 
numbers of M31 RGB stars are all located in the outer 
regions of M31's stellar halo (i?proj ^ 80 kpc), while the 
fields without kinematically-identified substructure span 
a similar range in projected distance from M31's center 
as the fields with kinematically-identified substructure. 
Fields with identified substructure are on average more 
metal-rich than fields without identified substructure, as 
would be expected if kinematically cold substructures are 
more metal-rich on average than the kinematically hot 
population (Fig. [21). 

3. SIMULATED STELLAR HALOS 

Simulations of hierarchical stellar halo formation pro- 
vide a physical explanation of the surface brightness ver- 
sus metallicity trend observed in our M31 halo data. In 
this section, we briefly describe the simulations (§[3T|), 
present the abundance properties of substructure in the 



simulated stellar halos fi) l3.2[) . and discuss the physical 
mechanisms behind the trends (§1331) ■ We will compare 
the abundance trends in the simulations with our M31 
observations in §|4| 

3.1. Methods 

The BJ05 simulations follow the formation of 11 MW- 
analog stellar halos via accretion in a cosmological con- 
text. For each host halo, a merger tree was gen- 
erated using the extended Press-Schechter f ormalism 
(|Somerville fc Kolattl fl999l : iLacev fc Coli [1993} . Only 
merger histories without recent high-mass mergers were 
selected, resulting in z = halos suitable for hosting a 
large disk galaxy like the Milky Way or M31. For each 
merger, an A'^-body simulation was used to track the evo- 
lution of the dark matter halo of the satellite galaxy in 
the host halo potential. The potential of the host galaxy 
is modeled with a smoothly growing, analytic function, 
and satellite-satellite interactions are neglected. Overall, 
these assumptions will lead to an unnatural enhancement 
of the level of substructure in the simulated halos, espe- 
cially in the inner halo. Thus, in the following analysis 
we discuss only the properties of the simulated stellar 
halos exterior to i? = 20 kpc. 

Gas accretion onto each satellite halo tracks its mass 
accretion history and the star formation rate within each 
object is proportional to its instantaneous gas content 
over a fixed timescale. The timescale is chosen so that 
dwarfs infalling today would contain similar gas fractions 
to those observed in Local Group field objects (see BJ05 
for details). Given the gas infall, star formation rates, 
a nd mass of a satel l ite as a function of time, the method 
of [Robertson et all (|2005f ) is used to follow t he chemical 
evolut ion of t he stellar populat i ons (s ee also [Font et all 
[2006bD . The 'Rob ertson et alj (|2005D prescriptions in- 
clude enrichment from Type la and Type II supernovae 
and the effects of feedback (from supernovae and stel- 
lar winds). The free parameters in the prescriptions are 
constrained by requiring consistency with observations: 
the Fe and a-element wind efficiencies are chosen sepa- 
rately by requiring matches to the observed stellar mass- 
metallicity re l ation for dwarf galax ies fe.g.. [Larsonlll974l : 
[Mated [19981 [Dekel fc Wool [2003[ ) and the a-element 
abundance patte rns in satellite dwarf spheroidals (e.g. 
[Venn et al.[ [20041 ) . Both gas accretion and star forma- 
tion are assumed to truncate once the satellite is accreted 
onto the host halo. A variable mass-to-light ratio is as- 
signed to every dark matter particle in order to embed 
the stellar components within the dark matter halos of 
the satellites in a way that reproduces the observed dis- 
tribution of structural properties of Local Group dwarfs. 
Each stellar particle in the simulations is randomly as- 
signed [Fe/H] and [a/F e] ratios based on the satellite's 
star formation history ([Font et al.[ [2006bl ). Thus each 
model satellite contains a realistic intrinsic metallicity 
spread, but no radial metallicity gradient. 

The resulting simulated halos have similar luminosities 
and stellar density profiles to the Milky Way's halo, and 
the surviving satellites are comparable in number, lumi- 
nosity distribution, and structural properties to those of 
the Milky Way (discovered prior to 2005). The mod- 
els are also able to reproduce the observed difference in 
[ck/Fe] ratios between the Milky Way's stellar halo and 
dwarf satellites. The lower [a/Fe] ratios in the surviving 
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simulated dwarf galaxies are due to their late accretion 
times and the prolonged, less bursty star formation his- 
tories of these system s compared with those of the pro- 
genito rs of the halo (jRobertson et al.l 120051 : iFont et all 
[20063) • In the models, the kinematically cold streams 
tend to have distinct chemical abundances (i.e. higher 
[Fe/H] and lower [g/Fel) than that of the smooth, more 
well- mixed halo ()Font et al.ll2006bl ). This is also due to 
the fact that the progenitors of these streams were ac- 
creted following the formation of the bulk of the halo. 

Of course, in the framework of the hierarchical forma- 
tion of stellar halos, there is no clear demarcation be- 
tween tidal streams and the "smooth" halo. In reality, 
present day stellar halos are expected to display a natural 
progession in the kinematical properties of their unbound 
substructure, ranging from thin and kinematically cold to 
more disperse, kinematically hot substructure. In light of 
the above findings, the prospect of using chemical abun- 
dances as a complementary dimension to phase-space to 
disentangle various substructures is encouraging. 

3.2. Predicted Abundance Properties of Substructure 

Figure [3] shows the projected surface brightness (top), 
metallicities (middle), and alpha-enrichment (bottom) for 
two of the BJ05 simulated stellar halos. The projected 
surface brightness reflects the total luminosity in each 
0.25 kpc X 0.25 kpc pixel, summed over all the stel- 
lar populations. The [Fe/H] and [a/Fe] values are the 
stellar-mass-wcighted average of the stellar populations 
in each pixel. These maps reveal that the highest sur- 
face brightness streams tend to be more metal-rich and 
less alpha-enhanced than the lower surface brightness 
streams. 

This is demonstrated quantitatively in Figure [5l which 
displays maximum surface brightness versus mean metal- 
licity (left panels) and alpha-enhancement (right pan- 
els) for individual stellar streams produced by accretion 
events in all 11 BJ05 s tellar halos, taken from the anal- 
ysis of I Johnston et all (2008. ) . Only stellar streams that 
contribute at least 50% of the total light in the halo at 
the position of their highest surface brightness point and 
which are at radial distances greater than 30 kpc from 
the center of the galaxy are included. Although for any 
given maximum surface brightness stellar streams exhibit 
a range of mean [Fe/H] and [a/Fe] values, the higher sur- 
face brightness streams are in general more metal-rich 
and less alpha-enhanced than the lower surface bright- 
ness streams. 

3.3. Physical Interpretation 

The points in Figure O are color-coded by mass of the 
progenitor satellite (top panels) and time of accretion of 
the progenitor satellite onto the host halo (bottom panels; 
see figure caption for details). The top left panel clearly 
shows the effect of the assumption of a mass-metallicity 
relation for the progenitor satellites (based on observa- 
tions of dwarf galaxies; ij |3.ip : more massive satellites 
are more metal-rich, thus the tidal streams produced as 
massive satellites fall into the host halo are also more 
metal-rich. Since the more massive progenitors also typ- 
ically produce brighter stellar streams, it follows that the 
surface brightness and metallicity of stellar streams are 
correlated. 
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Fig. 4. — Maps of projected surface brightness (top), aver- 
age metallicity (middle), and average alpha-enrichment (bottom) 
for 2 of the BJ05 simulated stellar halos. Each map depicts a 
300x300 kpc region. The [Fe/H] and [a/Fe] values are the stellar- 
mass-weighted average of all the stars in each pixel. A compar- 
ison of the surface brightness and metallicity maps reveals that 
the brightest features tend to be more metal-rich and less alpha- 
enhanced than the rest of the halo. Image credit: Sanjib Sharma. 

The alpha-enhancement of the stellar streams traces 
the time of accretion of the progenitor onto the host 
halo, as is demonstrated in the bottom right panel of 
Figure O Alpha elements are produced in Type II su- 
pernovae; their abundances become subsequently diluted 
by the elements produced in Type I supernovae. Thus, 
progenitors accreted early are alpha-enriched relative to 
progenitors that are accreted later and have had longer 
star formation histories. The surface brightness of stellar 
streams observed today is a function not only of the mass 
of the progenitor but also of the time since accretion, 
therefore the bottom right panel shows a slight trend for 
stellar streams produced by progenitors that were ac- 
creted at early times to have lower surface brightnesses 
than progenitors that were accreted at later times. This 
leads to the anti-correlation between surface brightness 
and alpha-enrichment. 

4. COMPARISON OF OBSERVATIONS AND 
SIMULATIONS 
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Fig. 5. — Maximum surface brightness versus average metallic- 
ity [left panels) and alpha enhancement {right panels) for stellar 
streams in the 11 BJ05 simulated stellar halos. The top pan- 
els are color coded according to the luminosity of the progeni- 
tor satellite [L < 10^ Lq (black); 10^ < L < lO'' Lq (red); 
10^ Lq < L < 10^ Lq (blue); 10* Lq < L < 10^ Lq (yellow)] while 
the bottom panels are color coded according to the time of accre- 
tion of the progenitor satellite onto the host halo [ 3 < tacc < 6 
Gyr (purple); 6 < iacc < 9 Gyr (cyan); 9 < tacc < 12 Gyr (or- 
ange)] . The distribution in the left panels is a reflection of the mass- 
metallicity relationship observed for Local Group dwarfs, while the 
distribution in the right panels reflects the time-sensitivity of alpha 
elements. The bottom panels also show a slight tendency for more 
recent accretion events to be brighter. 

4.1. Observing the Simulations 

Figure [5] provides useful insight into the physical mech- 
anisms underlying the expected general relationships be- 
tween the chemical enrichment and maximum surface 
brightness of stellar streams. The data in Figure [5] ap- 
pear to agree qualitatively with the characteristics seen 
in the left panels of Figure [51 identifiable (and thus 
fairly dominant and high surface brightness) substruc- 
tures tend to be relatively metal-rich. However, our 
Keck/DEIMOS spectroscopic masks provide information 
on the characteristics of the population in isolated lines- 
of-sight through M31's stellar halo, rather than global 
information about individual streams. The data cannot 
be compared directly with Figure [SJ since the maximum 
surface brightness of an observed stellar stream can only 
be determined by tracing the full extent of the stream in 
space. The only stellar stream for which this measure- 
ment could currently be made in M31 is the GSS. 

A more direct comparison between the data and sim- 
ulations can be made by analyzing stellar populations 
along individual lines-of-sight through the simulated stel- 
lar halos. Figure [6] displays mean stellar metallicity and 
alpha-enhancement as a function of surface brightness 
for sight lines through one of the BJ05 simulated stellar 
halos; each point represents a line-of-sight that covers a 
projected area of 0.4 kpc x 0.4 kpc. Only lines of sight 
beyond 20 kpc in projection from the center of the galaxy 
are included (§[SII])- Lines-of-sight through the simulated 
stellar halos whose stellar populations are dominated by 
debris from one or two accretion events are denoted by 
red solid circles in Figure [51 while lines-of-sight whose 
stellar populations are largely comprised of debris from 



-2.5 -2 -1.5 -1 -0 5 0.2 0.4 0.6 




-2.5 -2 -1.5 -1 -0.5 0.2 0.4 0.6 



[Fe/H] [a/Fe] 

Fig. 6. — Surface brightness versus stellar-mass- weighted aver- 
age metallicity (left) and alpha-enhancement (right) for individ- 
ual lines-of-sight through one of the 11 BJ05 simulated stellar ha- 
los, divided into two categories based on the number and promi- 
nence of individual stellar streams, (a) Lines-of-sight where the 
majority of the stellar population derives from one or two accre- 
tion events (solid red circles), (b) Lines-of-sight where the major- 
ity of the stellar population derives from multiple (> 3) accretion 
events (open blue squares) . ( c) The two types of sight-lines super- 
imposed. The overall trends of surface brightness with metallic- 
ity and alpha-enhancement (panel c), are similar to those in Fig- 
ureO although with a somewhat steeper slope. For a given surface 
brightness, lines-of-sight whose stellar populations are dominated 
by debris from one or two accretion events are on average more 
metal-rich and less alpha-enhanced than lines-of-sight whose stel- 
lar populations are composed primarily by debris from multiple 
accretion events. The surface brightness of the faintest tidal debris 
that has been identified as kinematically cold substructure in our 
Keck/DEIMOS survey is denoted by the dotted line (Fig. [2J while 
the current surface brightness limit of our spectroscopic survey is 
denoted by the dashed line (Fig.lSj. 

multiple (> 3) accretion events are denoted by blue open 
squares. 

The same general trends seen in individual streams 
(Fig. [5]) are repeated here, with a steeper slope: higher 
surface brightness lines-of-sight are on average more 
metal-rich and less alpha-enhanced than lines-of-sight of 
lower surface brightness. The same physical mechanisms 
discussed in ij |3.3l are relevant here, however each point in 
Figure [SI is averaging over stellar populations from multi- 
ple satellite progenitors. Thus, the trends in metallicity 
and alpha-enhancement with surface brightness are not 
as strong. 

4.2. Common Trends 

We can't know a priori the number of streams that are 
along any particular observed sight-line in M31's stellar 
halo. However, separating the kinematically cold and 
hot components in fields with identifiable substructure, 
as is done in Figure [2l effectively separates a single field 
into the equivalents of lines-of-sight dominated by stellar 
populations arising from a single accretion event (kine- 
matically cold components) and lines of sight dominated 
by stellar populations arising from many accretion events 
(kinematically hot components). In a stellar halo built in 
a purely dissipationless fashion, a population with a large 
spread in velocities is most likely to be formed from the 
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superposition of debris originating in multiple satellites 
that merged early with the host galaxy. The dotted line 
in Figure [S] denotes the surface brightness of the faintest 
tidal debris that is identifiable as kinematically cold sub- 
structure in our Keck/DEIMOS survey (Fig.[2|). 

The simulations (Fig. [6]) show the same general trends 
seen in our M31 data (Fig. [5]): (1) Higher surface bright- 
ness features are on average more metal-rich, and (2) 
at a given surface brightness lines-of-sight dominated by 
debris from one or two accretion events tend to be more 
metal-rich than those dominated by debris from multiple 
accretion events. 

We have discussed the first trend in §§ 13.21 - 1331 and 
shown that in a stellar halo formed via accretion, it 
can be explained via the mass-metallicity relationship 
for dwarf galaxies: larger galaxies are more metal-rich 
and also produce the brightest tidal streams. 

This same physical mechanism can also explain the sec- 
ond trend. In order for tidal debris from a single satel- 
lite to dominate the stellar population in a line-of-sight 
through a halo, the tidal debris must be part of a partic- 
ularly bright stellar stream relative to the average sur- 
face brightness at its location in the halo. It thus must 
originate from a relatively massive satellite (relatively 
metal-rich) and/or recent accretion event (low relative 
alpha-enhancement) compared to the accretion events in 
a line-of-sight of equivalent surface brightness but com- 
posed primarily of debris from multiple progenitors. 

With multiple physical mechanisms (e.g., the size and 
orbit of the progenitor and the time since its accretion 
onto the host halo) affecting the surface brightness of 
individual streams, a significant spread in the chemi- 
cal properties of individual sight-lines of a given surface 
brightness, as seen in Figure [6l is to be expected. This 
is particularly relevant for faint lines-of-sight, since faint 
stellar streams can be formed by the recent accretion of 
a low mass progenitor or early accretion events whose 
debris streams have spatially dissipated, or both. The 
brightest lines of sight, however, must be composed of 
debris from very massive and/or recent events, resulting 
in a narrower range in chemical properties. 

The data also show evidence of a narrower 
spread in metallicity among higher surface brightness 
streams/lines-of-sight and a larger spread in metallic- 
ity among fainter streams/lines-of-sight (Figs. [2] & [21 
i) l2.5|) . The rms deviation of mean [Fe/H] values for 
the high surface brightness {fiy < 29 mag arcsec"^) 
and low surface brightness (/iy > 29 mag arcsec"^) 
tidal streams observed in our M31 spectroscopic data 
are 0.08 dex and 0.32 dex, respectively (although the 
observed spread is likely underestimated due to system- 
atic errors, ij |2.5p . To compare the spread in [Fe/H] seen 
in the observations with that seen in the simulations, 
we impose surface brightness limits on the simulations 
based on the surface brightness of tidal streams observed 
in our survey (Fig. [2]), and measure the rms devation 
of the mean [Fe/H] values of lines-of-sight through the 
BJ05 simulated stellar halo which have a stellar popu- 
lation dominated by one or two accretion events (open 
circles in Fig. [6]). The rms deviation of [Fe/H] values 
for high surface brightness lines-of-sight (27 < fiy < 
29 mag arcsec"^) is 0.07 dex, while the rms deviation 
of [Fe/H] values for low surface brightness lines-of-sight 
(29 < hy < 31.7 mag arcsec^^) is 0.13 dex. The spread 



in [Fe/H] values in our observations and the simulations 
is similar for high surface brightness lines-of-sight, while 
the spread in [Fe/H] values in the observations is greater 
than in the simulations for lower surface brightness lines- 
of-sight. The rms deviation of [a/Fe] values in the simu- 
lations for the magnitude limits used above are 0.02 dex 
and 0.05 dex, respectively. 

We note that the absolute [Fe/H] scale is uncertain 
for both the observations (due to the age-metallicity de- 
generacy of the RGB as well as systematic errors in the 
generation of theoretical isochrones, ii l2.3[) and the simu- 
lations (where the output [Fe/H] is dependent on the in- 
put chemical enrichment prescriptions, § I3.1|) . Therefore, 
we do not attempt to directly compare the metallicities 
of the data and simulations. 

4.3. Discussion 

The common trend of increasing metallicity with 
increasing surface brightness in lines-of-sight through 
M31's stellar halo and simulated stellar halos is encour- 
aging. It gives us confidence that such simulations can 
lead to meaningful interpretation of stellar halo observa- 
tions. 

The surface brightness of M31's stellar halo is ob- 
served to decrease with in creasing distance fr om the cen- 
ter o f the galaxy ( Guhat hakurta et al.i 20()5l : llrwin et al.l 
I2005t llbata et al.ll2007l ). and the inner region of M31's 
stellar halo i s observed to be more metal- rich than the 
outer region (iKalirai et al.ll2006at IChapman et al.l [20061 : 
iKoch et all 1200 8') . ThusT^the general increase m mean 
[Fe/H] with increasing surface brightness of the kinemat- 
ically broad components in our M31 data (Fig. [21 open 
squares) is expected. In the hiearchical framework, this 
global trend in stellar halo metallicity could be explained 
with the same physical mechanisms invoked above fi^ l4.2p . 
Although the radial extent over which tidal debris from 
an individual accretion event is spread depends on the 
orbit of the progenitor, the effect of dynamical friction 
on the progenitor's orbit is more important for more mas- 
sive progenitors, causing them to sink to the center of the 
gravitational potential of the host galaxy more quickly 
than less massive progenitors. Therefore, more massive 
progenitors will tend to deposit most of their relatively 
metal-rich debris in the interior regions of the host halo. 

If the number of M31 stars detected in an observational 
field is small (< 20), it can be difficult to determine if the 
field is dominated by contributions from one or two kine- 
matical components or if the field contains contributions 
from multiple kinematical components (i.e., contains a 
relatively hot population). The results of the simula- 
tions can be used to examine the probable number of 
stellar streams in a given field. If the mean [Fe/H] of a 
field is relatively high compared to other fields of com- 
parable surface brightness, the field may be dominated 
by one or two stellar streams, while if the mean [Fe/H] 
is relatively low, the field is more likely to be dominated 
by multiple superposed stellar streams along the line of 
sight (i.e., the stellar population along the line of sight 
is approaching a dynamically hot population). 

There are several fields in our M31 spectroscopic sur- 
vey in which we have spectra of < 10 M31 RGB stars 
(open squares. Fig. [3]); all are in M31's sparse outer halo 
(^proj ^ 80 kpc). Based on the above discussion fi^ l4.2p . 
the outer halo fields with the highest mean [Fe/H] values 
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arc the most likely to be dominated by debris from one 
or two high mass progenitors, while the fields with lower 
mean [Fe/H] are more likely to be dominated by debris 
from many lower mass progenitors. Future observations 
will be able to test this conjecture. If fields with lower 
mean [Fe/H] are found to have stars with a wide range of 
velocities, it will support the hypothesis that the field is 
dominated by debris from many progenitors. Although 
the presence of a significant old, metal-poor stellar pop- 
ulation in a higher mass progenitor may somewhat ob- 
scure this effect, the mean [Fe/H] of a field dominated 
by such a progenitor would still be expected to be higher 
on average than the mean [Fe/H] of a field dominated by 
low mass progenitors which did not have an appreciable 
number of metal-rich stars. 

There are currently no observational constraints on the 
alpha-enrichment of stars in M31's stellar halo. How- 
ever, if a mean [a/Fe] can be measured for the M31 
RGB stars in a given field (i.e., by measuring the rel- 
ative abundances of elements in a composite spectrum 
formed by coaddition of spectra from multiple stars), it 
could be used to constrain the average time of accre- 
tion of the stellar population. Measurement of the mean 
[a/Fe] would greatly facilitate the interpretation of the 
surface brightness of specific fields by lessening the de- 
generacy between the effect of the luminosity of accreted 
satellites versus the time of accretion onto the host halo. 

5. SUMMARY 

We have presented the general trends of surface bright- 
ness with chemical enrichment of RGB stars observed in 
the course of our M31 Keck/DEIMOS spectroscopic sur- 
vey and of stellar streams in simulations of stellar halos 
formed in a hierarchical framework. Kinematically iden- 
tified substructure in M31's stellar halo is metal-rich rela- 
tive to the kinematically hot population, and higher sur- 
face brightness features are more metal-rich than lower 
surface brightness features. 

Similar trends are found in the simulated BJ05 stellar 
halos, which are formed entirely via accretion. Higher 
surface brightness lines-of-sight through the simulated 
stellar halos are in general more metal-rich than lower 
surface brightness lines-of-sight, and lines-of-sight dom- 
inated by stars from one or two accretion events are in 
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